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ABSTRACT 


In  recent  years  many  studies  in  fields  related  to  Meteor¬ 
ology  have  been  done  using  Oxygen-18  and  Deuterium  Isotopes.  These 
fields  include  Glaciology,  Paleoc limatology ,  Cloud  Physics,  and 
Hydrology.  In  this  thesis  some  of  the  local  processes  that  influence 
the  Oxygen-18  content  measured  in  precipitation  are  examined,  with 
the  ultimate  goal  of  using  Oxygen-18  as  a  tracer  of  the  history  of 
the  water  vapour  that  is  condensed  for  local  precipitation. 

Dansgaard  considered  the  mean  annual  del  values  in  precip¬ 
itation  of  coastal  stations,  and  found  that  they  were  linearly  re¬ 
lated  to  the  mean  annual  temperature  by  the  equation 

6  =  0.695  T  -  13.6  ppt 

By  applying  linear  regression  to  the  mean  monthly  del  values  and 
mean  monthly  temperatures  for  only  the  summer  months  (May  through 
September)  at  Edmonton,  Alberta  it  is  found  that  they  are  related  by 

6=0.  37  T  -  20.  0  ppt 

After  exploring  theoretically  some  of  the  physical  pro¬ 
cesses  that  cause  isotopic  fractionation  in  precipitation,  models 
are  developed  to  show  how  these  processes  modify  the  del  values 
measured  in  precipitation.  This  modelling  includes  both  saturated- 
pseudoadiabatic  cooling  and  evaporation  and  exchange  processes  acting 
on  the  raindrops  between  the  cloud  base  and  the  ground.  With  these 
models  it  is  shown  that  the  isotopic  fractionation  during  cooling  is 


■ 


dependent  upon  the  liquid  water  content  of  the  cloud  and  the  initial 
temperature  of  the  saturated  air.  This  fractionation  is  shown  to  be 
less  than  that  obtained  by  using  the  models  of  Miyake,  Matsubaya,  and 
Nishihara  who  considered  dry  adiabatic  cooling.  It  is  also  shown 
that  evaporation  and  exchange  can  cause  changes  of  the  order  of  two 
parts  per  thousand  under  summer  conditions. 

Grouping  the  data  by  individual  periods  of  precipitation 
showed  the  variations  in  del  values  of  precipitation  that  occur  with 
time.  It  is  found  that  different  types  of  storms  exhibit  different 
variations  of  del  values  with  time.  For  an  airmass  thunderstorm  the 
del  values  of  the  precipitation  are  found  to  decrease  markedly  during 
the  storm.  Attempts  to  explain  this  time  change  in  del  values  were 
made  using  the  thunderstorm  model  of  Byers  and  Braham  along  with  the 
models  that  were  formulated  in  this  study.  The  changes  in  time  of 
the  del  values  of  precipitation  during  general  periods  of  rain  are 
found  to  be  small  and  are  ascribed  to  evaporation  and  exchange  pro¬ 
cesses. 

In  this  study  it  was  found  that  the  data  available  were 
insufficient  to  permit  the  removal  of  the  local  changes  from  the 
measured  del  values.  It  is  recommended  that  further  sampling  be 
undertaken  so  that  future  work  may  establish  if  tracer  techniques 
are  feasible  using  Oxygen-18. 
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CHAPTER  1 


INTRODUCTION 

1 . 1  The  Isotopes  of  Oxygen 

In  nature  there  exist  three  stable  isotopes  of  the  oxygen 

atom.  They  have  atomic  weights  of  16,  17,  and  18  atomic  mass  units 

,  ,  ,  .  ,  16  17  18 

and  may  be  designated  as  0,  0,  and  0,  respectively.  The  most 

16  yi 

abundant  isotope  is  0  with  0,  having  an  odd  atomic  weighty  being 

the  least  abundant  in  natural  sources.  In  natural  waters  it  has  been 

1 8 

found  that  the  abundance  of  0  is  approximately  two  parts  per  thou¬ 
sand  (hereafter  designated  as  ppt)  (Dansgaard,  1953,  1964). 

1 . 2  Areas  of  Study  Related  to  Meteorology 

Oxygen  -18  concentrations  have  been  used  in  various  fields 
related  to  meteorology.  One  such  field  is  Glaciology.  Workers  such 
as  Gonfiantini  et  al.  (1963)  have  studied  the  variations  in  ^0  con¬ 
centration  that  occur  in  the  different  layers  of  a  glacier.  These 
variations  have  been  correlated  with  the  stratigraphy  and  have  shown 

annual  variations.  Epstein  and  Sharp  (1959)  used  the  variations  in 
18 

0  concentrations  with  position  in  the  glacier  along  with  a  model  of 

18 

glacier  movement  to  show  how  0  concentrations  could  be  used  in  glacier 
flow  studies. 

More  recent  studies  in  glaciology  by  Dansgaard  et  al.  (1969) 

and  Epstein,  Sharp,  and  Gow  (1970)  have  been  related  to  past  climate. 

These  workers  have  taken  deep  core  samples  through  the  ice  caps  of 

Greenland  and  Antarctica,  respectively.  Dansgaard  et  al.  (1969) 

1  8 

found  variations  in  0  concentration  dating  back  100,000  years.  In 
these  results  the  variations  due  to  the  Pleistocene  Ice  Age  and 


' 


2 

climatic  changes  since  that  time  were  shown.  Epstein  et  al.  (1970) 
found  similar  results  in  Antarctica. 

18 

Further  work  has  been  done  in  studying  the  0  concentration 
variations  in  hailstones.  Workers  such  as  Macklin,  Majoube,  and 
Merlivat  (1972),  Macklin,  Merlivat,  and  Stevenson  (1970),  and  Majzoub, 
Nief,  and  Roth  (1968)  used  a  model  to  explain  the  change  in  ^0  concen¬ 
tration  with  temperature  in  a  hail-producing  cloud.  With  this  model 

1 8 

and  the  observed  0  variations  in  the  hailstones  they  were  able  to 

infer  the  trajectories  of  the  hailstones  within  the  cloud. 

1 8 

Dansgaard  (1964)  studied  the  mean  annual  0  concentrations 
in  precipitation  at  stations,  for  which  data  were  available,  located 
around  the  North  Atlantic  Ocean.  By  comparing  these  mean  annual 
values  with  the  mean  annual  surface  temperature  for  the  station,  he 
developed  the  empirical  relation 


6  =  0.695  T  -  13.6  ppt 


(i.i) 


_  1  O 

where  6  is  the  annual  mean  of  a  measure  of  0  concentration  (to  be 
defined  in  Chapter  II),  and  T  is  the  mean  annual  temperature  in 
degrees  Celsius. 


1 . 3  Variation  in  Oxygen  -18  Abundance 

Oxygen  isotope  and  temperature  data  for  Edmonton,  Alberta 
were  compared  with  the  linear  relation  given  by  Eq.  (1.1).  Because 
annual  values  are  available  for  only  three  years  (1963  -  1965),  it  was 
assumed  that  mean  monthly  values  would  yield  the  same  relation.  Figure 
1  shows  the  mean  monthly  del  values  (6)  plotted  with  respect  to  the 
mean  monthly  temperature  data  collected  between  May  and  September  in 
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Mean  monthly  del  values  versus  mean  monthly  temperature 
,  Alberta  for  the  summer  months.  (X)  1961-1963,  (0)  1970- 
d  line  is  Eq.  (1.1),  dashed  line  is  Eq.  (1.2). 
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the  years  1961  to  1965  (from  International  Atomic  Energy  Agency,  1969, 
1970)  and  in  1970  -  1971.  It  is  immediately  apparent  that  all  of  the 
del  values  are  below  the  straight  line  value  of  Eq.  (1.1).  It  must  be 
noted,  however,  that  at  least  part  of  the  observed  difference  could  well 
be  due  to  the  change  from  considering  annual  values  to  considering 
monthly  values  for  only  the  summer  months.  Linear  regression  applied 
to  the  data  gave 

5  =  (  0.37  -  0.10  )  T  -  20  t  2.0  ppt  (1*2^ 


When  this  study  was  begun  it  was  thought  that  possibly 

1 8 

relations  such  as  Eq.  (1.1)  or  Eq.  (1.2)  could  be  applied  to  use  "  0 
concentrations  in  precipitation  as  a  tracer  to  study  the  source  of 
water  vapour  for  local  precipitation.  It  was  felt  necessary  to 

18 

consider  other  possible  mechanisms  to  explain  the  variations  in  0 
concentration  in  precipitation. 

Several  such  mechanisms  have  been  identified  in  the  litera¬ 
ture.  These  include: 

1)  the  isotopic  composition  of  the  original  source  of  the 
water  vapour, 

ii)  the  temperature  at  which  evaporation  and  condensation 
processes  occur, 

iii)  the  method  of  evaporation  from  the  source, 

iv)  the  past  condensation  history  of  the  vapour  between 
evaporation  at  the  source  and  condensation  into  precipitation, 

v)  the  effects  caused  by  the  mode  of  transport  (ie,  eddy 


diffusion  or  advection)  of  the  vapour  from  the  source  region,  and 


■ 
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vi)  exchange  and  evaporation  processes  acting  on  the  raindrop 
between  formation  and  analysis. 

Most  of  these  influences  will  be  examined  in  Chapters  III  to  V  after 

18 

considering  some  of  the  theory  of  °0  concentrations  in  Chapter  II,, 


CHAPTER  II 


THEORY 


2 . 1  Introduction 

18 

As  stated  in  Chapter  I  the  abundance  of  0  is  approximately 

two  ppt  in  natural  waters.  Absolute  concentrations  are  difficult  to 

1 8 

obtain.  However,  because  it  is  only  the  variations  in  0  content 

that  are  of  interest,  a  relative  change  in  content  is  measured.  Thus 

1 8 

the  relative  deviation,  or  del  value,  (denoted  by  £ )  of  the  0  content 
of  a  sample  from  a  reference  standard  is  measured.  Thus 

^  _  C  sample 

C  standard 

- 

where  C  is  the  ratio  of  mole  fractions  of  ^  -^0  to  H2  ^0  in  the  water 
(hereafter  referred  to  as  the  oxygen  ratio).  The  reference  standard 
used  is  SMOW  (Standard  Mean  Ocean  Water)  (Craig,  1961). 

If  a  sample  has  a  deviation  of relative  to  a  secondary 

standard  that  has  a  deviation  of  5^ from  SMOW,  Dansgaard  (1964)  gave 

s 

that  the  sample's  deviation  relative  to  SMOW  is 


x 


103  ppt 


(2.1) 


6  =  S'  +  S  +  (S'S  /103)  ppt 

s  s 

(2.2) 


2 . 2  Frac  tionation 

Isotopic  fractionation  of  water  occurs  whenever  water  under¬ 
goes  a  change  of  phase  from  vapour  to  liquid  or  solid  and  for  the 

reverse  change  of  phase.  This  is  due  to  H  ^0  having  a  higher  vapour 

1 8 

pressure  than  has  II ^  0.  A  fractionation  factor  a  is  defined  as 


a 


(2.3) 


. 
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where  C  and  C  are  the  oxygen  ratios  in  the  liquid  and  vapour  phases, 
i-j  V 

respectively.  If  the  vapour  pressures  of  H  and  are  denoted 

as  en ,  and  e  ,  respectively,  then  the  fractionation  factor  can  be 
lb  18 

approximated  by 


a  =  e  /e 

16  18 


(2.4) 


for  equilibrium  conditions.  Because  the  temperature  dependencies  of 

e1 ,  and  ein  are  not  the  same,  the  fractionation  factor  is  also  temp- 
16  18 

erature  dependent,  having  values  of  approximately  1.007  at  40  C  and 
1.014  at  -20  C  (Dansgaard,  1964). 


2. 3  Equilibrium  Processes 

To  demonstrate  how  the  oxygen  ratio  can  vary  in  space  and 

time  several  processes  will  now  be  considered. 

A  closed  system,  such  that  the  total  vapour  and  liquid 

phases  are  in  equilibrium,  will  now  be  considered.  If  a  small 

amount  of  vapour  is  condensed  in  this  closed  system,  the  newly- 

formed  condensate  has  an  oxygen  ratio  that  follows  Eq.  (2.3)  and 

1 8 

thus  the  vapour  becomes  depleted  in  0.  With  further  condensation 

both  the  vapour  and  the  newly-formed  condensate  become  more  depleted 

18„ 

in  0. 

Under  atmospheric  conditions  the  total  condensate  does  not 
stay  with  the  vapour,  but  falls  relative  to  the  vapour.  For  a  system 
that  is  closed  except  that  the  condensate  is  immediately  removed,  so 
that  no  exchange  occurs  between  the  liquid  and  vapour  phases,  Rayleigh 
condensation  occurs  (Dansgaard,  1964).  Then  for  the  condensation 


process 


. 


8 


(a  ~  1) 


dN 

v 

N 


v 


(2.5) 


where  N  is  the  mass  of  vapour  in  the  system.  As  shown  by  Dansgaard 

(1964)  this  leads  to  greater  fractionation  of  the  oxygen  isotopes, 

with  both  the  remaining  vapour  and  the  newly-formed  condensate  becoming 

1 8 

more  depleted  in  0  than  would  be  the  case  according  to  Eq.  (2.3). 

For  evaporation  of  a  fixed  finite  sample  of  liquid,  the  oxygen 
ratio  of  the  liquid  changes  according  to  an  equilibrium  batch  distil¬ 
lation  such  that 


(c.)  Fe 

L  o 


(2.6) 


where  (C  )  is  the  initial  oxygen  ratio  in  the  liquid,  F  is  the  fraction 
L  o 

of  the  liquid  remaining,  and 


a 


(2.7) 


(Craig,  Gordon,  and  Horibe,  1963,  Eq.  1).  As  F  ^  1  ande<0?Eq.  (2.6) 
leads  to  continued  enrichment  of  ^0  in  the  remaining  liquid  as  the  sample 
evaporates . 


2 . 4  Kinetic  Processes 

In  atmospheric  conditions  condensation  can  be  regarded  as 
a  slow  process  and  thus  equilibrium  conditions  as  discussed  above  will 
hold.  Evaporation,  however,  is  in  some  cases  a  fast  process  and  kinetic 
effects  must  be  considered.  These  influences  are  due  to  the  difference 
in  diffusivities  for  H2  ^0  and  H2  ^0  in  both  the  liquid  and  the  vapour 
phases.  These  differences  influence  the  fractionation  during  evaporation. 


. 
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Ehhalt  and  Knott  (1965)  derived  the  equation 


a  =  (t)l/2  (. _ 1 ±1 _ , 

K  n  In  ) 


'1-h'“DCV/CD 


(2.8) 


for  the  kinetic  fractionation  factor  o^^where  D**  and  D  are  the 
dif  fus ivit ies  of  Hg  ^0  and  H2  respect ively  in  air,  h'  is  the 
relative  humidity  of  the  air  measured  relative  to  the  drop  temperature, 
is  the  equilibrium  fractionation  factor  at  the  drop  temperature, 
and  Cy  and  are  the  oxygen  ratios  in  the  air  and  drop,  respectively. 


The  effective  fractionation  factor,  ct  is  then 

E 


aE  =  VaK 


(2.9) 


Because  a  is  usually  greater  than  unity,  the  use  of  the  effective 
K 

fractionation  factor  in  Eq.  (2.6)  leads  to  greater  distillation  than 
equilibrium  conditions. 


2 . 5  Equilibrium  Exchange  Effects 

A  droplet  of  water  in  air  containing  water  vapour  is  subject 
to  an  equilibrium  exchange  of  water  molecules  between  the  droplet  and 
the  vapour  even  when  no  net  mass  transfer  occurs.  Thus  the  droplet 
can  experience  a  change  in  the  oxygen  ratio  if  the  oxygen  ratio  of  the 
vapour  is  not  related  to  the  droplet's  oxygen  ratio  by  Eq.  (2.3)  when 
the  fractionation  factor  is  evaluated  at  the  droplet  temperature. 

For  this  exchange  Friedman,  Machta,  and  Soller  (1962)  and 
Miyake,  Matsubaya,  and  Nishihara  (1968)  developed  the  equation 


IIL 


m 


(nrn  -  nip)  exp  (-kt) 


(2.10) 
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where  m  is  the  amount  of  h/?  0  in  a  droplet,  subscripts  i  and  f  denoting 

initial  and  final  stages? respectively ,  t  is  the  duration  time  of  the 
exchange,  and  k  is  the  rate  constant  for  the  reaction  given  by 


3D0PM 

k  =  - 

r2ctpR*T 


(2.11) 


where  P  is  the  saturated  vapour  pressure  of  normal  water  at  temperature 

T  (in  this  and  further  equations  in  degrees  Kelvin),  R*  is  the  universal 

1 8 

gas  constant,  M  is  the  molecular  weight  of  H  0,  r  is  the  radius  of  the 
droplet,  p  is  the  density  of  the  droplet,  and  0  is  the  ventilation 
factor  of  a  falling  droplet  developed  by  Kinzer  and  Gunn  (1951), 

F  Re1/2 

0  =  1  + - -  (2.12) 

(4ttD/v)  !/2 


where  F  is  an  empirical  factor  of  Kinzer  and  Gunn,  Re  is  the  Reynolds 
Number,  and  v  is  the  kinematic  viscosity  of  the  air. 

Because  cloud  droplets  are  small,  giving  a  relatively  large 
rate  constant,  the  exchange  process  will  keep  the  cloud  droplets 
essentially  in  equilibrium  with  the  surrounding  vapour.  For  large 
raindrops  the  rate  constant  becomes  relatively  small,  so  that  longer 
adjustment  times  are  required  to  reach  equilibrium  with  the  surrounding 
vapour.  Thus  raindrops  are  seldom  in  equilibrium  with  the  vapout,  but 
they  are  modified  by  this  process. 


' 
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CHAPTER  III 


MODELLING 


3 . 1  Introduction 

In  an  attempt  to  explain  the  local  effects  on  the  oxygen 
ratio  that  appeared  in  the  data  collected,  a  computer  model  of  the 
theoretical  variations  of  the  oxygen  ratio  was  developed.  This  model 
included  both  adiabatic  condensation  of  water  vapour  and  modification 
of  the  oxygen  ratio  in  falling  raindrops  due  to  exchange  and  evaporation. 


3 . 2  Adiabatic  Condensation 

Miyake,  Matsubaya,  and  Nishihara  (1968)  developed  equations 

for  a  set  of  four  models  for  the  adiabatic  condensation  of  water  vapour. 

In  developing  these  equations,  however,  they  used  the  dry  adiabatic 

equation  for  the  change  of  volume  with  temperature,  and  thus  have 

neglected  the  latent  heat  of  condensation  released.  This  causes  large 

errors  in  the  amount  of  water  vapour  condensed  and  thus  in  the  change 

of  the  oxygen  ratio,  especially  at  higher  temperatures.  To  remove  this 

error  a  saturated-pseudoadiabat ic  cooling  process  was  used  to  obtain 

the  change  of  volume  with  temperature. 

Following  Miyake  et  al.  (1968)  four  models  were  developed. 

For  models  1  and  2  a  closed  system  was  considered,  so  that  all  condensed 

water  remained  in  the  system  during  further  cooling.  In  models  3  and  4 

an  open  system  was  considered  in  which  most  of  the  condensed  water  was 

removed  from  the  system,  except  that  the  cloud  retained  a  fixed  liquid 

water  density  W.  It  was  assumed  the  the  oxygen  ratio  was  a  close 

1 8 

approximation  to  the  0  concentration  (less  than  0.3  per  cent  error) 
and  thus  was  used  to  represent  both. 


■ 
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For  saturated  air,  the  mass  of  water  vapour,  N  ,  can  be 

obtained  by  the  expression  for  the  saturated  mixing  ratio,  w  , 

s 


Ntt  =  m  w 
V  s 


=  m 


ee 


P~e 


(3.1) 


where  m  is  the  mass  of  the  air  under  consideration,  e  is  the  ratio 

of  the  molecular  weights  of  water  and  dry  air  (=0,622),  e  is  the 

s 

saturated  vapour  pressure  with  respect  to  water,  and  p  is  the  total 
pressure.  To  obtain  the  saturated  vapour  pressure,  integration  of  the 
Clausius  -  Clapeyron  equation  for  evaporation  yields  . 

MT 

In  e  =  -  — — •  +  constant  (3.2) 

S  R~1 


where  M  is  the  molecular  weight  of  water,  L  is  the  latent  heat  of 
evaporation,  and  R*  is  the  universal  gas  constant.  Evaluating  Eq. 
(3.2)  at  0  C  (considering  L  to  be  constant  over  the  atmospheric  range 
of  temperature)  results  in 


e^  =  6.1078  exp  (A  +  ^) 

(3.3) 

where  e  is  now  the  saturated  vapour  pressure  with  respect  to  water, 
s 

in  millibars,  6. 1078  is  the  saturated  vapour  pressure  at  0  C  (List, 
1966)  and  the  constants  A  and  B  are  19.834  and  -5417.8,  respectively. 
For  saturated  air,  using  the  ideal  gas  law 

V  =  mRT*/p  (3.4) 


where  m  is  the  mass  of  air,  R  is  the  gas  constant  for  dry  air,  and  T* 
is  the  virtual  temperature  given  by 


13 


1  +  w  /  z 

T*  =  ( - - — )  T 

1  +  w 

s 


(3.5) 


The  fractionation  factor 
form 

a 


a  was  given  by  Miyake  et  al.  (1968)  in  the 


=  1  + 


(3.6) 


where  the  constants  c  and  d  have  values  of  -0.0105  and  5.8,  respectively. 


3.2.1  Model  1 

18 

For  a  closed  system,  the  mass  balance  equation  for  0  is 


(NvVo 


NVCV  +  nlcl 


and  for  the  total  water  substance  is 


(3.7) 


(Vo  =  NV  +  NL 


(3.8) 


where  (  )  denotes  the  initial  conditions,  N  is  the  mass  of  water 
substance,  C  is  the  oxygen  ratio,  and  subscripts  V  and  L  refer  to  the 
vapour  phase  and  liquid  phase,  respectively.  If  simple  equilibrium 
as  given  by  Eq.  (2.3)  is  considered,  then  Eqs.  (3.7)  and  (3.8)  can  be 
solved  to  give 


(Wo 


(-  -  1)  H  +  (Nv)o 


(3.9) 


1.  The  virtual  temperature  is  a  factor  used  in  the  ideal  gas  law  to 
correct  for  the  presence  of  water  vapour  in  the  air. 
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By  evaluating  and  a  according  to  Eqs .  (3d)  and  (3.6), respectively , 
the  oxygen  ratio  in  the  liquid  can  be  calculated  during  the  cooling 
process . 


3.2.2  Model  2 

Because  this  is  also  a  closed  system,  Eqs.  (3.7)  and  (3.8) 
remain  valid.  If  Rayleigh  condensation  is  considered,  Eq.  (2.5)  can 
be  integrated  to  give 


«Vi  =  (Vi-i 


<Vi  /  (Vi-i 


a-1 


(3.10) 


where  subscripts  i  and  i-1  refer  to  the  i  th  and  the  preceding  conden¬ 
sation  steps,  respectively.  Thus  can  be  determined  at  each  step. 
Combining  Eqs.  (3.7)  and  (3.8) 


^Wo  Vv 

(Vo  -  NV 


(3.11) 


and  the  oxygen  ratio  in  the  liquid  can  be  determined  during  the  cooling 
process . 


3.2.3  Model  3 

For  an  open  system  such  that  liquid  is  removed  from  the 
system,  the  conservation  of  mass  equations  (Eqs.  (3.7)  and  (3.8)  )  are 
no  longer  valid.  Thus  an  infinitesimal  condensation  was  considered. 
This  small  change  can  be  thought  of  as  a  closed  system,  and  the 
equations  for  the  conservation  of  the  oxygen  isotope  and  for  the 
water,  written  for  this  small  change,  become 
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d(NvCy)  +  d(NLCL)  =  0 


dN  +  dN 
V  Li 


=  o 


(3.12) 


If  simple  equilibrium  is  assumed  for  this  change,  substitution 
of  Eq.  (2.3)  and  Eq.  (3.12)  yields 


Ny  dCy  +  Cy  d  Ny  =  -  Nl  CCv  da  +  adCv)  +  aCv  d  Ny 


or,  on  rearrangement 


dC 


V 


'V 


(a-1)  dN  -  WVda 
N^  +  aWV 


(3.13) 


because  W  may  be  defined  as  N  /  V 


L 

Using  a  Taylor  Series  expansion  for  dN^T  and  neglecting  terms 


higher  than  second  order 


d(Nv). 

dT 


^Vi-f-1  ^Vi-1 

2  AT 


(3.14) 


Denoting  this  as  AN  ,  integration  of  Eq.  (3.13)  yields 


«Vi  =  (Vi-i  exp 


(a-1) AN  +  WVd 

Y7  (T.-T.  .) 
l  l  l-l 

Nv  +  WVa 


(3.15) 


that  can  be  evaluated  at  each  step.  By  using  Eq.  (2.3)  the  oxygen  ratio 
can  be  calculated  as  a  function  of  temperature. 


- 

16 


3.2.4  Model  4 

Assuming  Rayleigh  condensation  for  the  infinitesimal  change, 
Eq.  (2.5)  can  be  rearranged  to  yield 


+  cvdNv 


and  substituting  this  into  Eq.  (3.12)  yields 

aCy  dNv  =  -  Nl  i  q  +  CL  d  Nv 


or 


dC. 


C  -  aC 

—  N—  dNv 


(3.16) 


Because  Eqs .  (3.10)  and  (3.14)  apply 


dCL  CL  "  o(Vl-lt(Vi/(Vi-l1“’1  AN 


dT 


WV 


which  can  be  integrated  to  give 


(cl4 


[b+a  (C,).  ,] 

_ _ L  i-l  exp 


[a  (Ti-T±_1) ] 


b 

a 


(3.17) 


a  =  AN/WV 

where 

b  =  -aaCcy.^  [(Nv).  /  (Vt-d0"1 
Thus  the  oxygen  ratio  is  calculated  at  different  temperatures  during 
the  cooling. 

For  the  models,  the  pressures  were  obtained  from  a  tephigram 

for  a  saturated- pseudoadiabat  starting  at  the  lOOOmb  level  and  for 

every  -2  C  increment  in  temperature  along  the  curve  to  a  final  temperature 

of  -30  C.  From  Eqs.  (3.1)  and  (3.4)  the  amount  of  moisture  and  the 

volume,  respectively,  were  obtained  for  one  kilogram  of  dry  air  at  each 

1  8 

pressure- temperature  point.  These  data  along  with  an  initial  0 


. 


17 


concentration  in  the  vapour  supplied  the  basic  information  required  for 
the  four  models. 

3.2.5  Results 

To  demonstrate  the  results  of  the  saturated-pseudoadiabatic 
cooling  models  described  above,  water  vapour  with  an  oxygen  ratio  that 
was  in  equilibrium  with  SMOW  (a  del  value  of  0.0)  at  20  C  was  cooled  to 
-28  C  in  increments  of  2  C  for  each  of  the  models.  It  was  assumed  that 
super-cooled  liquid  water  was  the  condensed  phase  for  temperatures 
below  freezing. 

In  Figure  2  is  shown  the  variation  of  del  values  of  the  water 
with  cooling  for  models  1  and  2.  It  should  be  noted  that  although  the 
theory  of  Rayleigh  condensation  as  used  in  model  2  should  cause  a  greater 
fractionation  than  for  model  1,  the  imposition  of  the  closed  system  that 
retains  the  condensate  has  removed  the  effect  of  this  fractionation. 

As  has  been  mentioned,  a  closed  system  is  not  felt  to  be  a  realistic 
approach  to  the  atmosphere. 

Figures  3  and  4  show  the  variation  of  del  values  for  models 
3  and  4,  respectively,  assuming  liquid  water  contents  (denoted  as  LWC 
in  the  figures)  between  2.0  and  0.01  grams  per  cubic  meter.  These  open 
system  models  exhibit  much  larger  fractionations  than  do  the  closed 
system  models  discussed  above.  Close  examination  again  reveals  that  the 
equilibrium  model  3  has  a  greater  fractionation  than  does  the  Rayleigh 
condensation  model  4,  with  the  results  of  the  two  models  converging  as 
the  liquid  water  content  decreases  towards  zero.  This  is  again  due  to 
the  presence  of  the  liquid  water  in  the  system  that  allows  some  exchange 
and  thus  strictly  not  Rayleigh  conditions  that  require  the  immediate 
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TEMPERATURE  (C) 


Figure  2.  Del  values  as  a  function  of  temperature  for  pseudo- 
adiabatic  cooling  (models  1  and  2).  (b)  model  1,  (a)  model  2. 

Values  plotted  every  10  degrees  (calculated  every  2  degrees). 
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Figure  3.  Del  values  as  a  function  of  temperature  for  pseudo- 
adiabatic  cooling  (model  3).  LWC  of  2.0  (a),  1.0  (b) ,  0.5  (c) , 
and  0.01  ( d)  . 
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Figure  4.  Del  values  as  a  function  of  temperature  for  pseudo- 
adiabatic  cooling  (model  4).  LWC  of  2.0  (a),  1.0  (b),  0.5  (c), 
and  0.01  (d) . 


30 
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removal  of  the  condensate. 

Figure  5  demonstrates  the  changes  that  occur  if  different  in¬ 
itial  temperatures  are  used.  Initial  temperatures  of  30,  20,  10,  and  0 
C  were  used  with  the  vapour  in  equilibrium  with  SMOW  at  these  initial 
temperatures.  This  demonstrates  the  effect  of  the  latent  heat  released 
at  the  warmer  temperatures,  resulting  in  less  fractionation  for  the 
same  change  in  temperature  when  a  higher  initial  temperature  is  used. 

To  demonstrate  further  the  effect  of  the  latent  heat  released 
during  cooling  Figure  6  shows  a  comparison  of  results  of  the  models 
derived  above  (the  wet  models)  with  the  results  that  were  obtained 
using  the  models  derived  by  Miyake  et  al.  (1968)  (the  dry  models). 

The  same  initial  conditions  were  used  for  both  sets  of  models,  namely 
an  initial  temperature  of  30  C  and  water  vapour  with  an  content 
in  equilibrium  with  SMOW  at  that  temperature. 

For  models  1  and  2  the  wet  models  exhibited  del  values  that 
were  up  to  1.2  ppt  greater  than  those  given  by  the  dry  models,  a  change 
of  over  twenty-five  per  cent.  In  models  3  and  4  the  wet  models  give 
final  del  values  that  were  over  5.0  ppt  greater  than  those  given  by  the 
dry  models,  again  a  change  of  over  twenty-five  per  cent.  Other  liquid 
water  contents  for  models  3  and  4  showed  similar  results. 


3 „ 3  Evaporation 

For  an  evaporating  raindrop  Kinzer  and  Gunn  (1951)  developed 
the  equation 


=  4-irr  D(pa  -  p^)  9 


dM 

dt 


(3.18) 
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Figure  5.  The  effect  of  changing  the  initial  temperature  for 
model  3.  Initial  temperatures  of  0  C  (a),  10  C  (b) ,  20  C  (c) , 
and  30  C  (d).  LWC  of  1. 0  . 
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Figure  6a.  A  comparison  of  the  wet  models  to  the  dry  models  during 
cooling.  Solid  line  is  for  the  wet  models,  dashed  line  is  for  the 
dry  models.  Model  1. 
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Figure  6b. 


Model  2. 


» 
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Figure  6c. 


Model  3.  LWC  of  2.0 
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Figure  6d. 


Model  3.  LWC  of  0.5 
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Figure  6e.  Model  4.  LWC  of  0.5 
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where  M  is  the  mass  of  a  drop  of  radius  r,  p  is  the  vapour  density  at 
the  drop  surface  (subscript  a)  and  in  the  environmental  air  (subscript 
b),  0  is  the  ventilation  factor  defined  in  Eq.  (2.12)  and  D  is  the 
diffusivity  of  water  vapour  in  air,  given  by 


D  = 


(T) 


where  DQ  is  the  diffusivity  at  temperature  TQ  and  1000  mb  pressure 
(List,  1966). 


It  can  be  shown  that 


dM 

dt 


2-iTr  p 

w 


d(r2) 

dt 


where  P  is  the  density  of  the  water  drop,  and  thus  Eq.  (3.18)  becomes 
w 


d(r2) 

dt 


2D  (pa  '  V 


w 


(3.19) 


and  for  a  drop  falling  at  its  terminal  velocity 


Vr  (  =  dz/dt  ) 


d(r2) 

dz 


2D  (p  -  p  ) 
a  b 


V 

r 


P 


w 


(3.20) 


Substituting  Eq.  (2.12)  for  6  in  Eq.  (3.20)  and  expanding  Re  =  2rV  /v 


d  (r2) 


dz 


2D  (p  -  p  ) 
a  b 

V  p 
r  w 


1  +  F 


r  V 


2ttD 


1/2 


(3.21) 


Using  a  Taylor  Series  expansion  for  r(z  -  Az)  and  neglecting 


second  order  and  higher  terms 

2  ’  2 


r .  =  r 


d(r2) 


i-1  dz 


Az 


(3.22) 
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In  order  to  use  Eq.  (3.21)  it  was  assumed  that  the  temperature  of  the 
drop  was  equal  to  the  wetbulb  temperature  of  the  air.  The  terminal 
velocities  were  calculated  using  Eq.  (4)  of  Wobus,  Murray,  and  Koenig 
(1971).  This  equation  gives  an  empirical  fit  to  the  terminal  velocities 
determined  by  Gunn  and  Kinzer  (1949)  with  less  than  0.5  per  cent  error 
for  radii  larger  than  50  microns.  These  velocities  were  modified 
according  to  Eq.  (7)  of  Wobus  et  al.  (1971)  for  conditions  other  than 
a  pressure  of  1013.25  mb,  a  temperature  of  20  C,  and  a  relative  humidity 
of  50  per  cent.  Wobus  et  al.  (1971)  state  that  this  adjustment  gives 
terminal  velocities  up  to  10  per  cent  larger  than  those  measured. 

Using  the  results  of  Eq.  (3.22)  considering  different  environ¬ 
ments,  the  oxygen  ratio  of  the  evaporating  drop  was  calculated  using 
Eq.  (2.6).  Eqs.  (3.6)  and  (2.8)  were  used  to  calculate  the  fractiona¬ 
tion  factor  for  equilibrium  and  kinetic  effects,  respectively. 

3.  3.  1  Results 

To  demonstrate  the  effect  of  evaporation,  a  cloud  based  at 
1375  meters  above  the  ground  was  chosen.  The  surface  temperature  was 
26  C  with  a  dry  adiabatic  lapse  rate  to  the  cloud  base.  The  relative 
humidity  was  46.  6  per  cent  at-,  the  ground  and  increased  linearly  to 
100  per  cent  at  the  cloud  base.  The  water  vapour  below  cloud  base  had 
a  del  value  that  was  in  equilibrium  with  water  with  a  del  value  of 
-10. 0  ppt  at  20  C. 

Under  these  conditions  Figure  7  shows  how  the  del  values  of 
the  raindrops  changed  with  height,  based  on  an  initial  del  value  of 
-13.2  ppt  for  the  drops  and  with  equilibrium  evaporation  conditions. 

This  demonstrates  how  the  evaporation  effect  increases  with  decreasing 


u 

* 

■ 


- 


30 


6180 

Figure  7.  The  variation  of  del  values  of  raindrops  during  descent 
for  equilibrium  evaporation.  Initial  del  value  of  -13.2  .  Drop 
radii  of  1.0  mm  (a),  1.6  mm  (b) ,  and  2.5  mm  (c). 
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drop  radius  and  decreasing  relative  humidity,  with  a  total  change  of 
from  0.6  to  1.7  ppt  for  the  drops  considered. 

Figure  8  shows  that  the  effect  of  evaporation  on  the  del 
values  is  more  marked  if  kinetic  evaporation  is  considered,  even  for 
the  high  relative  humidities  experienced  near  the  cloud  base.  There 
were  total  changes  of  from  1.6  to  3.7  ppt  in  the  del  values  of  the 
drops  that  were  considered. 

3. 4  Exchange  Effects 

The  effect  of  isotopic  exchange  between  the  drop  and  the  en¬ 
vironment  on  the  oxygen  ratio  of  the  drop  was  calculated  using  Eq,  (2.  10) 

with  the  initial  drop  oxygen  ratio  as  nr  and  with  an  environmental 
liquid  concentration  as  given  by  Eq.  (2.3)  evaluated  at  the  drop  temp¬ 
erature  from  the  vapour  oxygen  ratio  as  m^.  The  elapsed  time  t  is 

given  by  t  =  z/V^  where  z  is  the  distance  fallen  at  velocity  Vr  . 

3.  4.  1  Results 

Equilibrium  exchange  was  considered  with  the  same  conditions 
as  described  in  section  3.3.1  .  Figure  9  shows  the  effect  of  exchange. 

It  is  noted  that  exchange  had  a  greater  effect  at  the  beginning  of 
the  process  than  near  the  end,  due  to  the  exponential  approach  to  an 
equilibrium  condition  with  the  vapour.  This  is  much  more  noticeable 
for  a  smaller  drop  (curve  a)  than  for  a  larger  drop  (curve  c) .  It 
should  be  noted  that  the  total  change  is  greater  than  that  for  equi¬ 
librium  evaporation  (Figure  7)  but  less  than  that  for  kinetic  evapor¬ 
ation  (Figure  8),  with  final  changes  of  from  1.3  to  2.6  ppt  for  the 
drops  considered.  Of  course  the  magnitude  of  the  effect  is  dependent 
upon  the  amount  of  disequilibrium  between  the  drops  and  the  vapour. 
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Figure  8.  The  variation  of  del  values  of  raindrops  during  descent 
for  kinetic  evaporation.  Lines  as  for  Figure  7. 
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Figure  9.  The  variation  of  del  values  of  raindrops  during  descent 
for  equilibrium  exchange.  Lines  as  for  Figure  7.  Vapour  in 
equilibrium  with  water  having  a  del  value  of  -10.0  at  20  C. 


CHAPTER  IV 


THE  DATA 


4. 1  Collection  Methods 

The  samples  used  for  this  study  were  collected  in  the  period 

from  June,  1970  to  September,  1971  in  the  vicinity  of  Edmonton, 
Alberta.  There  were  approximately  the  same  number  of  samples  for 
each  year. 

The  rain  samples  were  collected  in  open  pans  that  were 

placed  in  well  exposed  areas  during  periods  of  rain.  The  pan  was 

emptied  into  a  sample  bottle  that  was  then  closed.  Most  samples 

consisted  of  approximately  40  milliliters  except  when  the  rain  was 

too  brief  in  duration  or  too  light  in  intensity  to  yield  this  volume, 

so  that  smaller  samples  resulted.  In  some  instances  (during  very 

light  intensity  rain  or  when  using  a  small  collection  pan)  the 

sample  was  collected  over  a  long  period  of  time  in  an  attempt  to  get 

a  large  enough  sample.  In  these  cases  some  evaporation  from  the 

sample  could  have  occurred  before  it  was  collected.  This  evaporation 

1 8 

would  cause  an  enrichment  of  the  0  concentration  in  the  sample. 

Also  with  the  long  sampling  times  equilibrium  exchange  with  the  water 
vapour  in  the  air  could  have  modified  the  oxygen  ratio  of  the  sample 
if  it  differed  significantly  from  the  oxygen  ratio  of  the  vapour. 

The  snow  (hail)  samples  were  collected  from  areas  of  newly- 
fallen  snow  (hail)  and  were  placed  in  sealed  containers  before  being 
allowed  to  melt.  The  melt-water  was  then  transferred  to  a  sample 


. 
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bottle  to  be  kept  for  analysis. 

The  temperature  data  were  extracted  from  thermograph 
records  of  a  mesometeorological  network  in  the  city  of  Edmonton. 

The  temperatures  recorded  were  averaged  over  the  period  of  the 
sample  collection.  The  temperature  data  have  a  probable  error  of 
-1.0  G  considering  all  sources  of  error. 

4.  2  Analysis  Method 

The  samples  were  prepared  for  analysis  on  the  mass  spectro¬ 
meter  by  the  standard  method  as  given  by  Epstein  and  Mayeda  (1953). 

This  process  involves  equilibration  of  the  water  sample  with  carbon 

1 8 

dioxide  of  a  known  0  content  under  prescribed  conditions.  The 

2 

modified  carbon  dioxide  was  then  analysed  on  the  mass  spectrometer 

by  comparing  its  oxygen  ratio  to  that  of  the  standard  carbon  dioxide. 

The  preparation  of  the  sample  for  equilibration  was  carried 

out  in  a  very  low  pressure  (approximately  10  J  mm  of  Hg.)  system. 

Three  times  during  this  process  the  water  sample  was  frozen  to  very 

low  temperatures  (less  than  -140  C)  and  then  was  exposed  to  the  near 

vacuum.  If  on  any  of  these  occasions,  especially  the  first  when 

freezing  was  accomplished  with  the  sample  at  atmospheric  pressure, 

the  sample  was  not  completely  frozen,  evaporation  of  some  of  the 

sample  would  have  occurred.  This  evaporation  would  have  led  to  a 
1 8 

slight  0  enrichment  of  the  sample. 


Analyses  on  the  mass  spectrometer  were  carried  out  in  collabora¬ 
tion  with  the  Department  of  Physics,  University  of  Alberta. 


‘ 


- 


36 


Figure  10  shows  a  listing  of  the  computer  output  produced 

3 

during  analysis  of  one  sample.  The  numbers  2**7  and  2**2  specify 
that  128  consecutive  readings  were  averaged  to  produce  each  reading 
printed  and  that  an  average  of  4  readings  was  used  to  calculate  each 
trial  value,  respectively.  Four  trials  of  both  the  standard  (S)  and 
the  sample  (X)  were  run.  Then  the  ratio  of  the  sample  value  to  the 
standard  value  (X/S  ratio)  was  computed.  For  this  calculation  one 
of  two  alternate  methods  was  used.  If  a  central  sample  value  was 
considered  it  was  divided  by  the  mean  of  the  two  standard  values  that 
surrounded  it.  If  a  central  standard  value  was  considered  it  was 
divided  into  the  mean  of  the  two  sample  values  that  surrounded  it. 

The  two  methods  were  used  alternately  for  the  six  central  values  and 
thus  six  ratios  were  computed  and  printed,  along  with  the  mean  value 
and  the  standard  deviation  of  the  ratios.  If  the  standard  deviation 
was  greater  than  0.0002  (0.2  ppt),  one  or  two  X/S  ratios  could  be 
deleted  from  the  set  of  six  if  they  appeared  to  be  in  obvious  disa¬ 
greement  and  the  mean  and  standard  deviation  were  recalculated.  If  the 
standard  deviation  remained  too  large  another  three  trials  were 
performed  and  a  new  mean  X/S  ratio  and  standard  deviation  were  calc¬ 
ulated.  This  procedure  produced  an  X/S  ratio  (and  thus  del  value  by 
Eqs.  (2.1)  and  (2.2)  )  with  a  standard  deviation  of  0.2  ppt. 

This  measured  del  value  was,  strictly  speaking,  for  the 


J  Computer  program  was  written  by  R.  N.  Hunt  of  the  Department  of 
Physics. 
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INTTME=2**7 
GMT  I ME  =  2*  *2 
first  set  MUST  BE  s  ! 

THIS  SET  IS  S 

00*39788  00*39786  00*39779  00*39776 

00*39782  S 
THIS  SET  IS  X 

00*40006  00*40010  00*40009  00*40016 

00*40010  X 
THIS  SET  IS  S 

00*39810  00*39799  00*39802  00*39814 

00*39806  S 
THIS  SET  IS  X 

00*40017  00*40026  00*40019  00*40023 

00*40021  X 
THIS  SET  IS  S 

00*39812  00*39818  00*39815  00*39835 

00*39820  S 
THIS  SET  IS  X 

00*40032  00*40046  00.40047  00*40050 

00*40044  X 
THIS  SET  IS  S 

00*39843  00*39845  00*39848  00*39852 

00*39847  S 
THIS  SET  IS  X 

00*40072  00*40059  00*40063  00*40070 

00*40066  X 

01*00542  01*00526  01*00522  01*00533  01*00528  01*00521 

MEAN  RATIO  01*00528 
STD  DEV  IS  00*00008 
TYPE  Y  IF  DEV  SATISFACTORY?  Y 


Figure  10. 
sample. 


A  sample  of  the  computer  output  for  the  analysis  of  a 
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carbon  dioxide  that  was  equilibrated  with  the  samples.  Craig  (1957, 
pp  144-145)  gave  the  correction  factor  to  relate  this  measured  del 
value  to  the  actual  del  value  of  the  water  sample.  However,  because 
the  measured  del  values  were  small  and  the  correction  factor  is  near 
unity,  the  measured  del  value  was  assumed  to  be  that  of  the  water 
sample  (less  than  4  per  cent  error).  (See  Appendix  B  for  a  sample 
calculation  of  this  correction.) 

Table  A-l  in  Appendix  A  gives  a  listing  of  all  the  samples 
that  were  collected.  This  includes  the  del  values  of  the  sample,  the 
air  temperature,  the  del  value  of  the  standard  used,  the  X/S  ratio 
measured,  and  the  location  of  the  collection  site. 


CHAPTER  V 


DISCUSSION  AND  RESULTS 


5 . 1  Introduction 

In  discussing  the  results  of  analysis  of  the  precipitation 
data  collected,  only  the  summer  period  (the  months  of  May  through  Sept¬ 
ember)  will  be  considered.  Furthermore,  only  those  samples  for  which 
surface  air  temperature  observations  were  available  were  included. 

The  data  were  grouped  by  month  and  monthly  means  of  temperature 
and  del  value  were  calculated.  Table  1  lists  the  number  of  samples  in 
each  month,  the  mean  surface  air  temperature,  and  the  mean  del  value 
and  standard  deviation  for  each  month. 

Figure  11  shows  these  means  plotted  and  includes  the  straight 
line  of  Eq.  (1.1)  from  Dansgaard  (1964).  As  in  Figure  1  all  mean  del 
values  fall  below  those  given  by  Eq.  (1.1)  as  would  be  expected. 

5 . 2  General  Results 

In  Chapter  I  several  causes  for  a  variation  in  del  value  in 
precipitation  were  given.  Some  implications  from  these  will  now  be 
developed. 

As  can  be  seen  in  Figures  1  and  11  and  in  Table  1  the  mean 
del  values  for  precipitation  in  the  area  of  interest  were  near  -16  ppt 
for  the  summer  months.  Thus  if  this  water  was  the  source  of  water 
vapour  for  the  precipitation  and  assuming  equilibrium  evaporation  and 
condensation  processes,  the  first  condensate  would  have  had  a  del  value 
of  near  -16  ppt  and  further  condensation  would  have  yielded  lower  del 
values  and  the  raindrops  in  the  cloud  would  have  had  a  del  value  slightly 
lower  than  -16  ppt.  As  shown  in  Chapter  III  exchange  and  evaporation 


' 
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Figure  11.  Mean  monthly  del  values  versus  mean  temperature  for  the 
collected  samples.  Horizontal  bars  indicate  one  standard  deviation 
about  the  mean.  Line  is  Eq.  (1.1).  Samples  marked  1 T '  are  summer  means. 
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TABLE  1 


A  summary  of  the  samples  collected  during  the  summer  months 


\  YEAR 

1970 

1971 

\ 

MONTH 

NO 

T 

6 

<r 

NO 

T 

6 

<T 

ppt 

ppt 

ppt 

ppt 

MAY 

0 

- 

- 

- 

6 

13.3 

-13.4 

0.9 

JUNE 

3 

8.9 

-18.3 

2.1 

44 

12.0 

-19.1 

1.6 

JULY 

29 

16.0 

-14.4 

1.6 

13 

15.2 

-16.9 

1.8 

AUGUST 

9 

12.5 

-17.1 

2.0 

3 

17.2 

-15.2 

1.1 

SEPTEMBER 

7 

13.2 

-11.0 

1.5 

4 

8.9 

-18.5 

1 . 6 

TOTAL 

48 

14.5 

-14.7 

1.8 

70 

12.7 

-18.0 

1.8 

. 
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would  have  raised  this  value  so  that  the  precipitation  would  have  had 
a  del  value  of  near  -16  ppt.  Also  as  will  be  shown  below  (section  5„4) 
evapotranspiration  raises  the  initial  del  value  of  the  vapour  and  thus 
of  the  precipitation.  This  local  evaporation  is  most  marked  for 
convective  precipitation. 

Another  possible  large  source  of  the  water  vapour  for  precip¬ 
itation  was  the  Pacific  Ocean.  It  may  be  considered  to  have  a  del  value 
of  near  0  ppt.  Thus  for  equilibrium  evaporation  and  condensation  the 
first  condensate  would  have  had  a  del  value  of  near  zero.  Thus  the  rain¬ 
drops  would  have  higher  del  values  than  those  considered  above.  However, 
to  get  to  the  area  of  interest  this  moisture  must  have  crossed  the  mountain 
barrier  between  the  ocean  and  Alberta.  During  this  traverse  of  the 
mountains  some  of  the  moisture  will  have  precipitated  out  and  thus  only 

part  of  the  original  water  vapour  comes  across  the  barrier. 

1  8 

As  0  is  preferentially  condensed  out  due  to  the  fraction¬ 
ation,  this  precipitation  on  the  mountains  preferentially  removed  part 
1  8 

of  the  0.  Thus  when  the  remaining  water  vapour  reached  Alberta  it  was 
18 

depleted  in  0  with  respect  to  its  original  content  and  thus  further 

precipitation  in  Alberta  would  have  had  del  values  such  as  those  measured. 

18 

This  preferential  loss  of  x  0  through  precipitation  may  be  regarded  as 
a  reservoir  effect.  Further  changes  during  movement  across  the  mountains 
can  be  caused  by  the  mode  of  transport  as  described  in  detail  by  Eriksson 
(1965). 

5 . 3  The  Change  in  Del  Values  with  Time 

An  alternate  method  of  examining  the  data  was  to  consider 
the  change  in  del  values  with  time  during  a  period  of  precipitation. 
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In  the  data  collected  there  were  ten  different  sets  that  included  at  least 
five  sequential  samples.  Each  of  these  sets  was  then  plotted. 

5.3.1  Instability  Precipitation 

Figure  12  shows  the  change  of  del  values  with  time  that 
occurred  for  an  airmass  type  thunderstorm.  It  was  noted  that  del 
values  decreased  with  time,  with  a  change  of  near  five  ppt.  There 
are  several  possible  explanations  for  this  effect. 

One  possibility  is  the  effect  of  exchange  and  evaporation. 

At  the  beginning  of  the  storm  the  air  below  the  storm  would  have  been 
much  drier  and  the  cloud  base  higher  than  near  the  end  of  the  storm. 

Thus  evaporation  of  the  raindrops  would  have  been  much  more  marked 
at  the  initial  stages  and  raised  the  del  values  of  the  initial  time 
and  thus  the  del  values  would  have  decreased  with  time.  The  effect 
of  exchange  would  be  similar. 

A  second  possibility  is  the  reservoir  effect  that  was 
decribed  in  the  previous  section  (5.2).  This  would  have  caused  a 
general  decrease  of  the  del  values  with  time  as  the  cloud  became 
depleted  in  water  vapour.  A  further  explanation,  the  one  that  the 
author  believes  was  the  most  probable,  will  be  explained  in  detail 
in  the  following  section  (5.4)  . 

In  Figure  13  is  shown  the  change  in  del  values  that  occurred 
with  a  thunderstorm  associated  with  an  upper  level  disturbance  that 
released  the  instability,  rather  than  surface  heating  having  released 
the  instability .  It  is  immediately  noted  that  there  was  no  general 
decrease  with  time  of  the  del  values  as  was  noticed  in  the  airmass 
thunderstorm.  If  it  is  assumed  that  little  surface  air  was  drawn  into 
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Figure  12.  Time  variations  in  del  values  of  precipitation  associated 
with  an  airmass  thunderstorm  on  July  23,  1971.  Vertical  barbs  on 
lines  show  one  standard  deviation  for  measured  del  values. 
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Figure  13.  Time  variations  in  del  values  of  precipitation  associated 
with  a  non-airmass  thunderstorm  on  August  13,  1970. 
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the  thunderstorm  the  effect  described  in  Section  5.4  would  have  been 
absent.  The  higher  del  values  at  the  beginning  and  the  end  of  the 
storm  were  due  possibly  to  the  evaporation  effect  as  was  described 
previously. 

5.3.2  General  Precipitation 

Figure  14  shows  the  change  in  del  values  with  time  that 
occurred  in  a  synoptic  situation  that  produced  warm  frontal  rain. 

As  the  front  did  not  move  appreciably  towards  the  area  of  interest, 
the  general  increase  in  del  values  attributable  to  a  type  of  reservoir 
effect  (Dansgaard,  1953)  did  not  take  place.  Rather  the  decrease 
that  was  evident  at  the  beginning  of  the  observations  can  be  attributed 
to  evaporation  effects  that  decrease  with  time  as  the  cooler  airmass 
beneath  the  warm  air  became  more  and  more  saturated. 

Figures  15  or  16  show  the  results  for  cold-low  type  of 
precipitation.  Figure  15  shows  nearly  constant  del  values  with  time 
as  there  existed  no  transitory  effects  to  cause  the  del  values  to  change 
The  almost  continual  rain  situation  shown  in  Figure  16  also  shows  near 
constant  del  values  over  a  period  of  days.  The  higher  del  value  during 
the  daylight  hours  on  the  15th  was  probably  due  to  an  evaporation  effect 
both  of  the  drops  as  they  fell  and  of  the  partial  sample  in  the  collect¬ 
ion  pan  because  the  precipitation  at  this  time  was  very  light  in  inten¬ 
sity  and  intermittent  in  time.  Exchange  effects  between  the  sample  in 
the  collection  pan  and  the  water  vapour  in  the  air  could  have  given  the 
nearly  constant  del  values  due  to  the  long  sampling  times.  The  cause 
of  the  higher  del  values  at  the  end  is  undetermined. 
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Figure  14.  Time  variations  in  del  values  of  precipitation  associated 
with  warm  frontal  rain  on  September  1,  1970. 


48 


Figure  15.  Time  variations  in  del  values  of  precipitation  associated 
with  cold  low  precipitation  on  June  5,  1971. 
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HOUR  OF  DAY 


Figure  16.  Time  variations  in  del  values  of  precipitation  associated 
with  cold-low  precipitation  on  June  14-18,  1971. 
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5 . 4  An  alternate  Approach 

It  has  been  shown  that  in  shower-type  precipitation  the 
oxygen  ratio  decreases  with  time.  As  mentioned,  part  of  this  decrease 
with  time  can  be  attributed  to  the  reservoir  effect  in  the  cloud. 
Additionally,  the  cloud  conditions  can  be  used  to  provide  a  further 
insight  into  the  observed  variations. 

Byers  and  Braham  (1948)  proposed  a  model  for  the  life 
cycle  of  a  heavy  shower  or  thunderstorm  cloud.  The  model  consisted  of 
three  stages;  the  cumulus  stage,  the  mature  stage,  and  the  dissipating 
stage.  In  the  cumulus  stage  the  cloud  exhibited  updrafts  in  all 
regions  with  cloud  droplet  growth  but  with  no  precipitation.  In 
the  mature  stage  there  existed  a  zone  of  downdrafts  associated  with 
heavy  precipitation,  with  updrafts  maintained  in  other  areas  in  the 
cloud.  The  final  dissipative  stage  of  the  cloud  contained  down- 
drafts  throughout  with  lighter  intensity  precipitation  falling. 

In  order  to  model  oxygen  isotope  variations  in  precipitation 
from  a  cloud  of  this  type  an  atmosphere  having  a  constant  wetbulb 
potential  temperature  of  20  C  was  assumed.  It  was  furthur  assumed 
that  in  the  mature  stage  the  zone  of  downdrafts  that  contained  the 
heavy  precipitation  could  be  represented  by  saturated-pseudoadiabatic 
cooling  with  a  liquid  water  content  of  two  grams  per  cubic  meter  and 
a  cloud  base  at  approximately  fourteen  hundred  meters  (cloud  base 
temperature  of  12  C) .  Similarly  the  dissipative  stage  was  repre¬ 
sented  by  saturated-pseudoadiabatic  cooling  with  a  liquid  water 
content  of  one-half  gram  per  cubic  meter  and  a  cloud  base  at  approx¬ 
imately  five  hundred  meters  (cloud  base  temperature  of  16  C) .  The 


. 
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del  value  of  the  cloud  droplets  at  0  C  in  the  cloud  was  presumed  to 
be  representative  of  the  del  values  of  the  raindrops  when  they  reached 
the  base  of  the  cloud  for  both  stages. 

In  all  calculations  an  initial  vapour  del  value  of  -19.  1 
ppt  was  considered.  This  corresponded  to  the  vapour  del  value  if 
water  with  a  del  value  of  -10.0  ppt  was  evaporated  under  equilibrium 
conditions  at  20  C,  or  alternatively  to  water  with  a  del  value  of 
-19.1  ppt  that  underwent  evapotranspi ration.  Eriksson  (1965)  stated 
that  during  evapotranspiration  no  fractionation  can  occur  in  this 
process  because  the  very  minute  size  of  the  capillary  network  of  the 
evaporating  medium  permits  no  mixing  to  occur. 

The  effect  of  changing  the  liquid  water  content  of  the 
cloud  proposed  for  the  mature  stage  is  shown  in  Figure  17.  The 
lines  drawn  are  for  liquid  water  contents  of  2.0,  1.0,  and  0.5  grams 
per  cubic  meter.  A  change  of  0.8  ppt  at  0  C  was  noted.  The  variation 
caused  by  changing  the  cloud  base  temperature  for  a  constant  liquid 
water  content  is  shown  in  Figure  18.  At  0  C  a  change  of  2.6  ppt 
occurred  for  the  conditions  shown  in  this  example.  Figure  19  shows 
the  results  of  combining  these  effects  of  liquid  water  content  and 
cloud  base  temperature.  This  yielded  a  total  change  of  -2.3  ppt 
between  the  conditions  chosen  for  the  mature  stage  and  for  the  dissi¬ 
pating  stage. 

1 8 

Further  modification  of  the  0  content  of  the  precipitation 
due  to  equilibrium  exchange  and  evaporation  between  the  cloud  base 
and  the  ground  must  also  be  considered.  Figure  20  shows  the  final 
del  values  for  drops  of  different  radii  for  the  atmospheric  conditions 
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Figure  17.  The  effect  of  liquid  water  content  on  the  del  values  of 
cloud  droplets  during  cooling.  Model  3.  LWC  of  2.0  (a),  1.0  (b) ,  and 
0.5(c).  Cloud  base  temperature  of  12  C. 
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Figure  18.  The  effect  of  cloud  base  temperature  on  the  del  values  of 
cloud  droplets  during  cooling.  Model  3.  Cloud  base  temperatures  of 
12  C  (a),  16  C  (b),  and  20  C  (c).  LWC  of  2. 0  . 
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Figure  19.  The  combined  effect  of  liquid  water  content  and  cloud 
base  temperature  on  the  del  values  of  cloud  droplets  during  cooling. 

T  of  12  C,  LWC  of  2.0  (a);  T  of  16  C,  LWC  of  1.0  (b) ;  and  T  of  20  C, 
LWC  of  0.5  (c). 


55 


Figure  20.  Final  del  values  of  raindrops  due  to  equilibrium  ex¬ 
change.  Solid  line  for  mature  stage,  dashed  line  for  dissipating 
stage. 
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described  above,  examining  only  the  equilibrium  exchange  with  the 
vapour  in  the  air  below  cloud  base.  This  resulted  in  a  final  del 
value  difference  of  near  2.7  ppt  between  drops  of  the  same  radius  in 
the  mature  stage  compared  to  drops  in  the  dissipative  stage. 

Figure  21  displays  the  results  of  equilibrium  evaporation 
for  the  raindrops.  Because  the  mature  stage  has  a  higher  cloud  base 
and  a  drier  atmosphere  below,  the  drops  experience  a  larger  evapor¬ 
ation  than  for  the  dissipating  stage.  This  effect  is  more  marked 
for  the  smaller  drops  as  a  larger  percentage  of  these  drops  is 
evaporated.  The  final  del  values  show  differences  of  between  2.8 
and  3.  8  ppt  for  the  radii  considered.  Kinetic  evaporation  is  shown 
in  Figure  22.  This  process  yielded  even  larger  differences,  ranging 
up  to  4.5  ppt.  A  summary  of  the  results  is  shown  in  Table  2.  This 
table  includes  the  effect  of  changing  the  initial  del  values  of  the 
drops.  The  values  shown  are  the  change  in  the  final  del  value  expres¬ 
sed  as  a  percentage  of  the  change  in  the  initial  value.  These  values 
are  dependent  on  the  atmospheric  conditions  below  the  cloud  base,  and 
are  shown  for  both  the  mature  stage  and  the  dissipating  stage. 

Figure  23  shows  this  change  in  percentage  with  radius  for  both  stages. 

The  modelling  as  described  in  Chapter  III  did  not  allow  for 
combining  these  results  of  exchange  and  evaporation.  A  conservative 
estimate  of  combining  them  gives  a  difference  of  3  to  4  ppt,  as  the 
effects  of  evaporation  will  slightly  reduce  the  effects  of  exchange, 
but  not  vice-versa.  Thus  this  proposal  is  seen  to  produce  results 
of  the  right  order  of  magnitude. 


. 


. 


57 


TABLE  2 

Summary  of  the  changes  of  del  values  that  occurred  below  cloud  base, 
including  the  effect  of  a  change  of  initial  drop  concentrations. 
Values  shown  are  the  final  del  values  of  the  drops,  and  the  change 
that  occurred  between  the  two  stages. 


PROCESS 

FINAL  DROP  DEL  VALUE 

BELOW  CLOUD  BASE 

| EFFECT  OF  CHANGING  INITIAL 
CONCENTRATION 

| 

MATURE  D IS  S I PAT ING 

RADIUS  STAGE  STAGE  CHANGE 

MATURE  DISSIPATING 

STAGE  STAGE 

mm  pp  t 

|  7c  7c 

EQUILIBRIUM 

EXCHANGE 

1.0 

-10.5 

-13.1 

-2.6 

21.7 

57.8 

1.6 

-11.3 

-14.1 

-2.7 

45.6 

74.8 

2.5 

-11.8 

-14.5 

-2.7 

59.5 

82.4 

EQUILIBRIUM 

EVAPORATION 

1.0 

-11.5 

-15.3 

-3.8 

100 

100 

1.6 

-12.3 

-15.4 

-3.1 

100 

100 

2 .5 

-12.6 

-15.4 

-2.8 

100 

100 

KINETIC 

EVAPORATION 

1.0 

-9.5 

-14.0 

-4.5 

61.8 

79.4 

1.6 

-10.9 

-14.6 

-3.7 

76.1 

87.9 

2.5 

-11.5 

-14.9 

-3.4 

82.8 

91.7 

' 
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Figure  21.  Final  del  values  of  raindrops  due  to  equilibrium  evapor¬ 
ation.  Lines  as  for  Figure  20. 
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Figure  22.  Final  del  values  of  raindrops  due  to  kinetic  evaporation. 
Lines  as  for  Figure  20. 
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Figure  23.  Percentage  of  the  initial  change  in  del  value  of  rain¬ 
drops  exhibited  by  the  final  raindrops.  Solid  lines  are  for 
exchange,  dashed  lines  for  kinetic  evaporation.  (X)  mature  stage, 
(0)  dissipating  stage. 


chapter  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


6. 1  Conclusions 

In  Chapter  I  it  was  shown  that  mean  monthly  data  for 

Edmonton,  Alberta  did  not  fit  Dansgaard's  (1964)  regression  (Eq. 

(1.1)  ),  but  rather  followed  Eq.  (1.2).  It  must  be  again  noted  that 

Eq.  (1.1)  was  obtained  using  mean  annual  values  while  Eq.  (1.2) 

was  obtained  using  mean  monthly  values  for  only  the  summer  months. 

Again  assuming  that  similar  regressions  would  result  using  either 

method,  two  conclusions  can  be  drawn.  Firstly,  the  difference  in 

constants  between  the  two  relations  implies  that  precipitation  in 

1 8 

Alberta  is  depleted  in  0  with  respect  to  the  coastal  stations  used 
in  determining  Eq.  (1.1).  This  depletion  may  be  regarded  as  a 
continentali ty  effect  (caused  by  a  type  of  reservoir  effect). 
Secondly,  the  decreased  slope  of  Eq.  (1.2)  compared  with  Eq.  (1.1) 
can  be  attributed  to  a  change  in  the  type  of  precipitation.  Thus 
while  coastal  stations  experience  mainly  precipitation  associated 
with  organised  weather  systems,  continental  stations  such  as  Edmonton 
receive  more  instability  type  precipitation  during  the  summer  months 
and  thus  experience  a  wider  range  of  types  of  precipitation  than  do 
the  coastal  stations. 

It  has  been  shown  that  by  making  certain  simplifying 
assumptions  it  is  possible  to  explain  part  of  the  variations  that 
occurred  in  the  data  collected.  However,  because  of  the  assumptions 
made,  it  was  possible  to  show  only  the  qualitative  effects  and 
general  trends  for  the  causes  of  variations  listed  in  Chapter  I. 
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Thus  it  is  not  possible  to  remove  the  variations  due  to  local  effects 

1 8 

from  the  data  to  enable  the  0  to  be  used  as  a  tracer  for  the  water 
vapour. 

Studies  of  changes  in  del  values  with  time  within  individual 
precipitation  periods  have  shown  several  features.  In  most  cases 

18 

studied  the  effect  of  evaporation  enriching  the  precipitation  in  0 

for  at  least  part  of  the  time  has  been  demonstrated.  Also  demonstrated 

to  a  lesser  degree  was  the  effect  of  exchange.  By  the  use  of  a  cloud 

1 8 

model  the  marked  decrease  in  time  of  0  concentration  associated  with 
airmass  thunderstorms  was  demonstrated. 

6. 2  Recommendations  for  Further  Studies 

The  nearly  constant  del  values  with  time  measured  in 
general  precipitation  do  not  exclude  the  possibility  of  exchange 
during  collection  affecting  the  samples.  Thus  the  collection  methods 
must  be  reviewed  and  improved  to  minimize  the  effects  of  evaporation 
during  collection  and  of  exchange  with  the  water  vapour  in  the  air. 

In  this  way  a  check  of  the  observed  constancy  in  the  precipitation 
del  values  can  be  made  to  see  if  this  is  a  natural  phenomenon. 

The  practical  use  of  the  models  developed  requires  some 
method  of  sampling  the  water  vapour  in  the  air.  Thus  some  method 
of  collecting  these  samples  must  be  developed.  The  collection  of 
cloud  droplet  samples  would  help  to  complete  the  observational  net¬ 
work,  although  this  could  prove  to  be  a  monumental  task.  It  is  also 

recommended  that  Hydrogen  isotopic  analysis  of  the  samples  be  under- 

18 

taken,  because  the  relation  between  Deuterium  and  0  concentrations 
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gives  a  measure  of  the  relative  importance  of  the  kinetic  processes 

that  occur  (Dansgaard,  1964). 

The  models  that  have  been  developed  may  prove  useful  in 

other  fields  of  study.  As  mentioned  in  Chapter  I,  some  work  has 

1 8 

been  done  in  the  last  few  years  using  0  in  the  study  of  hailstones 

with  models  similar  to  the  closed  models  developed.  It  is  felt  that 

the  use  of  the  open  models  developed  would  be  more  accurate  for  this 

work.  By  expanding  and  further  developing  the  ideas  expressed  for 

18  . 

convective  storms  it  may  prove  possible  to  use  0  as  an  independent 
check  for  numerical  models  for  clouds  that  are  being  developed  for 
other  uses. 

1 8 

While  at  the  present  time  it  is  impossible  to  use  0  as 
a  tracer  for  the  water  vapour,  it  is  hoped  that  the  models  developed 
can  be  used  in  future  studies  with  a  more  complete  set  of  data  to 
investigate  this  possibility. 
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Figure  A-l.  Time  variations  in  del  values  of  precipitation 
associated  with  an  airmass  thunderstorm  on  July  20,  1970. 
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Figure  A-2.  Time  variations  in  del  values  of  precipitation 
associated  with  an  airmass  thunderstorm  oh  July  20,  1970. 
(Different  observing  site  than  for  Figure  A-l. 
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Figure  A-3.  Time  variations  in  del  values  of  precipitation 
associated  with  a  non-airmass  thunderstorm  on  July  9,1970. 
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Figure  A-4.  Time  variations  in  del  values  of  precipitation 
associated  with  a  period  of  general  rain  on  June  14-16,  1971. 
Dashed  and  solid  lines  indicate  different  observing  sites. 
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Figure  A-5.  Time  variations  in  del  values  of  precipitation 
associated  with  a  period  of  general  rain  on  July  29,  1970. 
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TABLE  A-2 
The  site  code 


CODE  LOCATION 

01  11504  -  3 7 B  AVENUE 

02  11030  -  74  AVENUE 

03  112  ST.  *  SASKATCHEWAN  DRIVE 

04  11130-57  AVENUE 

05  1 7 A  MEADOWLARK  VILLAGE 

06  4020  -  112  AVENUE 

07  SYLVAN  LAKE 

08  POPLAR  BAY 

09  PIGEON  LAKE 

10  SULPHUR  MOUNTAIN 

99  UNKNOWN 


TABLE  A-3 


The  weather  (WX)  code 


CODE  WEATHER 

00  WATER  SAMPLE 

60  LIGHT  INTERMITTENT  RAIN 

61  LIGHT  CONTINUOUS  RAIN 

6  2  MODERATE  I  NT E  PM  I TT  ANT  RAIN 

63  MODERATE  CONTINUOUS  RAIN 

70  LIGHT  INTERMITTANT  SNOW 

71  LIGHT  CONTINUOUS  SNOW 

80  LIGHT  RAIN  SHOWER 

81  MODERATE  OR  HEAVY  RAIN  SHOWER 

89  HAIL  SHOWER 

95  THUNDERSHOWER 

96  THUNDERSHOWER  WITH  HAIL 

97  HEAVY  THUNDERSHOWER 

99  HEAVY  THUNDERSHOWER  WITH  HAIL 


The  equilibration  of  the  carbon  dioxide  and  water  during 
sample  preparation  causes  a  change  in  the  oxygen  ratio  of  both  the 
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carbon  dioxide  and  the  water  due  to  the  finite  masses  used.  To 
correct  for  this  modification  Craig  (1957,  p  144)  gave 

6=5  IP±gl  -  —  6  (B.  1) 

r  m  p  p  s 

where  is  the  corrected  del  value  relative  to  the  standard  with  a 
del  value  of  5  ,  is  the  measured  del  value  during  analysis,  and 

a  is  the  fractionation  factor  for  the  equilibration  process,  given  by 


a 


C  / 

carbon  dioxide 


water 


(B.  2) 


where  C  is  the  oxygen  ratio  in  the.  carbon  dioxide  and  in  the  water,  and 

has  a  value  of  1.039  at  25  C.  p  is  the  ratio  of  gram  atoms  of  oxygen 

in  the  water  sample  to  gram  atoms  of  oxygen  in  the  carbon  dioxide 

introduced  into  the  sample  flask. 

During  preparation  of  the  samples  1  to  2  ml  of  water  were 

placed  in  a  25  ml  flask  that  was  then  filled  with  carbon  dioxide  to 

a  pressure  of  approximately  38  cm  of  mercury.  Thus  p  has  values  of 

between  55.56  and  115.2;  the  values  for  1  and  2  ml  of  water,  respectively. 

With  6  having  a  value  of  -19.3  ppt,  Eq.  (B. 1)  becomes 
s 

6  =  1.019  5  +  0.361  ppt 

r  m 

for  1  ml  of  water,  and 

<5  =  1.009  6  +  0.174  ppt 

r  m  ^ 


. 


for  2  ml  of  water.  Considering  an  average  of  1.5  ml  of  water  (  p 
equals  84.88),  Eq.  (B. 1)  becomes 
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5r  1.012  6^  +  0.236  ppt  (B. 3) 

This  correction  means  that  all  del  values  given  must  be 

increased  by  0o 24  ppt  and  that  the  actual  del  values  have  a  larger 

spread  than  the  measured  del  values. 

Using  one  of  the  larger  measured  del  values  (see  Table  A-l, 

sample  no.  JL  005)  6  equals  7.  18,  and  assuming  a  small  volume  of 

m 

water  used  (1  ml),  6  equals  7.68,  giving  a  maximum  error  of  7  per 
cent  in  the  relative  del  value,  or  an  error  of  4  per  cent  in  the  del 
value  with  respect  to  SMOW  (-11.8  versus  -12.3),  Using  Eq.  (B.3) 
this  becomes  a  likely  error  of  4  per  cent  (  6  equals  7.43)  in  the 
relative  del  value,  or  a  likely  error  of  2  per  cent  in  the  del  value 
with  respect  to  SMOW  (-12.0  versus  -12.3). 

Because  most  samples  had  a  smaller  measured  del  value,  the 
error  introduced  by  neglecting  this  correction  factor  was  smaller. 
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